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Abstract. We summarize the main results of a broad analysis on electrostatic, spherically 
symmetric (ESS) solutions of a class of non-linear electrodynamics models minimally coupled 
to gravitation. Such models are defined as arbitrary functions of the two quadratic field 
^ ■ invariants, constrained by several physical admissibility requirements, and split into different 

families according to the behaviour of these lagrangian density functions in vacuum and on the 
boundary of their domains of definition. Depending on these behaviours the flat-space energy 
{3J^' of the ESS field can be finite or divergent. For each model we qualitatively study the structure 

of its associated gravitational configurations, which can be asymptotically Schwarzschild-like or 
with an anomalous non Schwarzschild-like behaviour at r — )■ oo (but being asymptotically flat 
and well behaved anyhow). The extension of these results to the non-abelian case is also briefly 
Q.^ [ considered. 

o 
o 



1. Introduction 

' The non-linear Born-Infeld (BI) model [Ij, introduced in 1934 as an attempt to classically 

provide a finite energy for the electron field, has been the cornerstone of many developments in 
theoretical physics in the last three decades. The BI lagrangian density takes the form 

LBi=p'\i-\ii-i,-i^A, (1) 




where X = -iF^j.F^'^ = - and Y = -\F^^F*^''' = 2E ■ H are the two quadratic field 
invariants that can be constructed with the field strength tensor F^jy = d^Aj^ — d^A^ and its 
dual F*^^'^ = Fa0, and where E and H are the electric and magnetic fields, respectively. 

When minimally coupled to gravity, its electrostatic spherically symmetric (ESS) solutions 
have been studied for a long time [2]. With no doubt, the finding that string theory leads to 
(abelian and non-abelian) Bl-like lagrangians in its low-energy limit [3], has greatly contributed 
to the renewed interest in this kind of non-linear actions. As a result many non-linear 
electrodynamics (NED) models, aside from BI one, leading to asymptotically flat ESS structures, 
have been considered in the literature in the last decades Moreover, the characterization of 
the full class of physically admissible NED theories leading to asymptotically Schwarzschild- 
like configurations when coupled to gravity was performed in Ref.[5]. Other models have 
been found to lead to well behaved gravitational configurations but approaching asymptotic 



flatness slower than the Schwarzschild term [6l [?]. These lines of research have also lead to 
some promising attempts in the obtention of regular (singularity-free) ESS black hole solutions 
in NEDs minimally coupled to gravity [8]. Let us finally mention that black hole solutions 
generated by NEDs in other gravitational backgrounds have been also considered j9]. 

In this work we summarize the main results of a general procedure for the classification and 
proper characterization of asymptotically fiat electrically charged black hole solutions of NEDs 
defined as arbitrary functions of the two quadratic field invariants without explicitly fixing the 
lagrangian density function. The key point in such an analysis is the existence of a first-integral 
of the field equations, which takes the same form with and without gravity. These models 
are constrained by a number of physical admissibility requirements and classified according to 
the asymptotic and central-field behaviours of their ESS solutions. Their associated fiat-space 
energy can be finite or divergent, due to one or both of these field behaviours in the limits of the 
integral of energy (r — )• oo and r ~ 0). For each one of these models we qualitatively characterize 
their associated (asymptotically fiat and well behaved) gravitational configurations. Finally, the 
extension of these procedures to the non-abelian case is briefiy discussed. 



2. Gravitating ESS fields 

Let us begin by introducing the action for NED models minimally coupled to the Einstein action 
as 

S = [ d^x^a f - ^{X, Y)) , (2) 



VievrG 

where ip{X, Y) is an arbitrary function of the two field invariants previously introduced. We 
consider an ansantz for the static spherically symmetric metric as 

ds2 = e^(^)(it2-e'^(^)(ir2-rW, (3) 

where = dO^ + sit? 6d(fP' is the interval on the 2-sphere. However, due to the symmetry of 
the NED energy-momentum tensor T/ = T^, only a single independent metric function remains 
as A(r) = e'^^'") = e~^(^). The Einstein equations G^^, = SttT/^^ (we are taking units G = 1) for 
this metric can be easily obtained to give [5] 

A(0 = l-^ + ^^, (4) 

where M is an integration constant, identified as the Arnowitt-Deser-Misner (ADM) mass, and 
we have introduced, by convenience, the exterior integral of energy, defined as 

eeAr,q) = / R^T^{R,q)dR. (5) 

However, the procedure leading to Eq.dH) needs the implicit assumption that the NED source is 
such that eex{f,Q) converges in its upper limit (r — )• oo). As we shall see in the sequel there are 
NED models not satisfying this property, and in such cases the metric is integrated as 

A(r) = l + ^-^^^, (6) 
r r 

where now C is an integration constant not longer identified as the ADM mass (since, as we 
shall see, the large-r behaviour is governed in this case by the electrostatic field term, rather 
than the "Schwarzschild" one) and now 

e^n{r,q) = r R^T^{R,q)dR (7) 
Jo 



is the interior integral of energy. In this case, the obtention of Eq.([6]) needs the assumption that 
^inif^Q) converges in the lower hmit (r — >■ 0). But even if neither s^x 

(r, g) nor ein{r,q) can be 

defined, then the Einstein equations can be stiU integrated as 



>/ ^^ D 2e(r,q,0) 

A(r,g,Z)) = l + --^^;^, (8) 



eir,q,B)=47r r''Tl^{r,q) + B, (9) 



where 



B and D = C — 2B being integration constants. Consequently, in order to properly characterize 
the metric function in the different cases, one is lead to analyze the conditions under which 
^ex{i'-,q) and ein{r,q) are well defined. This amounts to study the behaviour of E{r) at r — )• oo 
and around r ~ 0. To do this we take into account the fact that the NED field equations 
(we are defining ipx = ||), obtained as V^{ipxFi"' + (pyF*'"') = lead, for ESS fields 
[E{r) = E{r)^ / 0, H{r) = 0), to a first-integral 

rVx^(r) = q, (10) 

where q is an integration constant identified as the electric charge associated to the field E{r). 
Remarkably, this first-integral takes the same form, in the Schwarzschild coordinate system, as in 
the absence of gravity in spherical coordinates. Consequently, we can perform the classification 
of the models in fiat space according to the asymptotic and central-field behaviours of E{r), and 
Eq. (fTO|) will allow the translation of the results obtained there to the gravitational context. 



2.1. Definition of the models 

Let us now establish a set of constraints on the NED models, leading to what shall be called 
admissible models. We restrict the lagrangian density function (p{X, y = 0) to be continuous, 
derivable and single-branched in its domain of definition of the X — Y plane, which must be open, 
connected and including the vacuum (X = Y = 0). The condition (p(X,Y) = ip{X,—Y) must 
be also satisfied in order to implement parity invariance. In addition we require the positive 
definiteness of the flat-space energy density 

P = T^ = 2X^x - v>{x) > (Vx2 + y2 + ll + y 1^ - ^{x, y) > o, (ii) 

whose fulfillment implies the conditions v?(0, 0) = 0, ^|(x>o,y=o) > and Lp{X > 0,0) > (see 
[lOj). These admissibility conditions endorse the definiteness and single- valued character of the 
ESS solutions and are essential in the analysis of the problem considered here. In particular, 
from the first-integral (jlOp it follows that E{r) must be a monotonic function of r for a given 
q. Moreover, for admissible models £ex{f^q) is a monotonically decreasing and concave 
function of r, while Sin{r,q) is a monotonically increasing and convex one. 
For finite-energy ESS solutions, the energy associated to the ESS field is given by 

e{q) = / r^T^{r,q)dr = q^/^e{q = 1), (12) 
Jo 

where e{q = 1) is the energy of the solution of unit charge, a universal constant for a given 
model. We can now split the families of admissible models into two sets, according to the 
character (finite or divergent) of the energy ()12p . Such a character only depends, for admissible 
models, on the behaviour of the ESS fields as r — )• oo and at r = 0. In this sense, we shall 
consider power-law expressions for the field behaviour around each of these limits. Let us then 
study each case separately. 



I) Models with finite-energy ESS solutions 



Asymptotically the field must behave as 



E{r) ~ rP, (13) 

with p < —1 for convergence of the energy there. The behaviour of the associated lagrangian 
density is obtained, by making use of the first-integral (fTO|l . as 

ip{X)r^X^. (14) 

In what follows we shall call this asymptotic behaviour, given by Eq. (jl3p . as B-type fields. Such 
behaviour can be slower than Coulombian (—2 < p < —1, subcase Bl), Coulombian (p = —2, 
B2) or faster than Coulombian (p < —2, B3). 

As r ~ there are two ESS field behaviours compatible with the finite energy requirement 
(A-type fields) for admissible models. In case Al the field behaves at the center as in Eq. ijlSp 
but now with — 1 < —p < for the convergence of the energy there (and, consequently, the field 
diverges in this case), and the lagrangian behaves also as in Eq. ()14p . On the other hand, in case 
A2 the field becomes finite at the origin, behaving there as 

E{r) ^ a- br", (15) 

where a and cj(> 0) are fixed parameters for a given model while 6 is a function of the charge 
[5]. Concerning the lagrangian density, its behaviour around the central values of the fields is 
given by (the case a = 2 is singular, see Ref.[10]) 

(^(X,y = 0)~^^^^(a-^/X)^+A, (16) 
2 — a 

where A = (p{a'^,Y = 0) is an integration constant. Consequently we see that there are two 
different behaviours for (p{X): if o" < 2 the lagrangian density diverges at X = a'^ while if o" > 2 
it takes a finite value there given by A. 

Consequently, for any admissible model (p{X, 1" = 0) to support finite-energy ESS solutions, 
these fields must belong to one of the B-cases at r — )• oo, and to one of the A-cases around 
r ~ 0. In between both domains the lagrangian density must be a strict monotonically increasing 
function of X. 

II) Models with divergent-energy ESS solutions 

ESS fields whose energy diverges at r = are obtained from Eq. ()13p . but now imposing the 
restriction p < —1 there. We shall call this kind of behaviour as UVD (meaning "ultraviolet" - 
divergent). A similar behaviour is assumed in the asymptotic limit, with — 1 < p < for 
(asymptotically vanishing) ESS solutions with divergent energy there. The acronym IRD 
("infrared" -divergent) will be employed for calling this family. Consequently the class of 
admissible NED models with divergent-energy ESS solutions is naturally classified into three 
families, resulting from the combination between i) divergent-energy fields at the center but not 
asymptotically (UVD-B), ii) asymptotic divergent-energy but not at the center (A-IRD) and 
iii) divergent energy in both limits (UVD-IRD). As a consequence of these behaviours Eexif, q) 
can only be defined in the first one of these cases, and em(r, q) only in the second one. 

2.2. Metric structure. Horizons 

With these ingredients we resume the previous discussion on the metric structure. First we look 
for the horizons r^ of the different configurations, which are the solutions of A(r) = in Eq.([3]) 
(in such cases where e^x 

{r,q) can be defined, see the previous section), implying 



(p(X, Y = 0) 







A2 f ( 
(CT<2)/, 


Al 
-l<p<0) 


-2<p<-l 








/^2 ~r/ 

(a>2yi/ 




/p = -2 






B3 






p<-2 






B2 

/ 












Z^^^:^^^^ IRD 















X 



Figure 1. Qualitative behaviour of the famihes of admissible lagrangian density functions 
^p{X^Y = 0). The behaviour at X ~ (r — )• oo) is split into finite-energy cases (Bl, B2 
and B3) and the energy-divergent one (IRD). At the center of the ESS solutions, there are 
two A2 cases (X ~ a^): for cr < 2, ^piX) diverges there while if o" > 2, (/7(-'^) becomes finite. 
Al is a vertical parabolic branch with divergent slope for large X. UVD is the family of 
divergent-energy solutions there, with its different possible slopes plotted. 



M-^ = eex(r,g), (17) 

and thus horizons are given by the cutting points between the monotonically decreasing and 
concave curve Sexij"^ q) the beam of straight lines M — |, corresponding to different values of 
the ADM mass M (see Fig. 2). Let us now study, following this procedure, each class of models, 
resulting from the combination of one of the possible field behaviours at the center (UVD, Al 
or A2) with one of the asymptotic cases (Bl, B2, B3 or IRD). 
I) Al-B : At r — 7- oo the metric (jl]) vanishes as {p < —1) 

2M 327rg „ , , 

A r) ^1 + -j^ -rP + ■■■ (18 

^ ' r (p-2)(p + l) ^ ^ 

and, consequently, the large-r behaviour is governed by the Schwarzschild term, regardless of the 
B-subcase chosen (Bl, B2 or B3). At the center eex{f,q) takes the value e((?), with a negative 
divergent slope there. Consequently, the cutting points with the beam M — | give the following 
configurations (see Fig. 2): i) There is a family of extreme black holes (EBH), satisfying the 
conditions A(r) = ^^^^^1^^ = 0, which give the EBH radius and mass as 

8T^rl^xtrToirhextr,q) = 1 ; Mextr{q) = ^"^y^ + Ao{rhextr, Q) , (19) 

where Ao{r,q) is the electric potential; ii) M < Mhi^xtr^Q)'- naked singularity (NS); iii) 
MfiextriQ) < M < £{q): two-horizons (Cauchy and event) BH; iv) M > £{q)- single-horizon 
BH; v) M = £{q): for this critical configuration we must study the behaviour of the metric 
around r ~ ( — 1 < p < 0): 



(20) 



Figure 2. Qualitative behaviour of eex(^i l) for the famihes Al, A2 (three subcases) and UVD 
at the center combined with B-fields at r — )• oo, cut by a set of dashed hues, corresponding to 
different values of M. The cut points give the horizons of the different configurations (see the 
discussion in the text). All the finite-energy cases have been chosen in order to begin at the 
same value of e{q) as r — )• 0. 



and we see that in this case A(r) diverges there to — oo, as in subcase iv) and, consequently, 
this is a single-horizon BH as well. In subcases iv) and v) the central singularity is spacelike, 
as opposed to the timelike singularity of subcases i), ii), and iii). Example of this family: The 
Euler-Heisenberg lagrangian ip{X^Y) = "f- + l^i^"^ + |^^)) ^ being a constant (see Ref.[5] for 
details of this family in gravitation) . 

II) A2-B : Now there are three subcases, according to the (now finite) slope of eexi'>',Q) 
around r ~ 0. A2a (IQirqa > 1): Here there are five classes of configurations: i) M < Mhextrio)'- 
NS; u) M = Mhextriq)- EBH; iii) Mhextr{q) < M < e{q): two-horizons BH; iv) M > e{q): single- 
horizon BH; v) M = e((/): for this configuration the metric at the center is finite, behaving there 
as (fj > 0) 

2(M-e(g)) 327rbq ^ . o 

r (cr +l)(2-(7) 

and, consequently, in this case A(0) = 1 — IGirqa < 0, leading to a single-horizon BH as well. 
A2b (IGvrga < 1): in this branch there are not EBH, and the available configurations are i) 
M < e{q): NS; ii) M > e{q): single-horizon BH; in) M = e{q): in this case A(0) = l-lQirqa > 0, 
as can be seen from (j21|) . leading to a NS. A2c (167rga = 1): the structures now are similar as 
in case Al-B, excepting for M = £{q), since A(0) = and this configuration is a kind of "black 
point" [5j. Example of this family: The BI lagrangian defined by Eq.([T]). 

III) UVD-B : In this case Eex{T~,q) diverges as r — )• (since now £{q) becomes divergent, 
as explained in subsection 2.1.) and there are three possible configurations: i) M < MhextriQ)- 
NS; ii) M = Mh(,xtr{<l)'- EBH; iii) M > Mh(.xtr{Q)'- two-horizons BH. Example of this family: 
The Maxwell lagrangian ip{X,Y = 0) = aX (a a constant), leading to the Reissner-Nordstrom 

solution A(r) = 1 - ^ + ^. 

IV) Al (and A2)-IRD: As already discussed, in this case Eexif^o) cannot be defined and 
the metric is integrated as in Eq.©. Now horizons are determined by the equation 



-^=e^n{r,q), (22) 

and by taking into account the monotonically increasing and convex clictrBjCter of (r, q) the 
different configurations can be classified fohowing the same procedure as in the previous cases, 
now using the cutting points with the beam of straight hues ^^^y^ in Eq. ()22p . corresponding to 
different values of the integration constant C [7]. For the Al and A2a cases at the center we 
find: i) A family of EBHs, whose radius is given by that of Eq. ()19p . now with an associated 
value C = Cextriq)', ii) C > Cextrio)- NS; iii) < C < Cextrio)- two-horizons BH; iv): C < 0: 
single-horizon BH; v) C = 0: single-horizon BH (now the metric diverges to — oo in case Al, 
while takes a finite value A(0) = 1 — IGvrga < in case A2a). In cases A2b and A2c the 
extreme black hole and two-horizons BH configurations disappear and we are lead to i) C > 0: 
NS; ii) C < 0: single- horizon BH; iii) C = a NS in case A2b and an extreme black point [7] 
in case A2c. This class of gravitational configurations approaches asymptotic flatness slower 
than the Schwarzschild solution, since the asymptotic behaviour is governed by the ESS term 
(see Eq. p8|) ). Example of this family: The lagrangian (p{X,Y = 0) = X'^ with 3/2 < 7 < oo 
(Al family at the center). 

V) UVD-IRD : In this case neither eexii",q) nor ein{r,q) exist, but the above procedure 
for classifying the different gravitational structures works again, employing now Eqs-dH]) and 
Q, together with the value of the integration constant D. This leads to: i) as in the 
previous cases, EBH configurations are found for a value D = DhextriQ)] h) D > Dh^xtriQ)- 
NS; iii) D < DhextriQ)'- two- horizons BH. Example of this family: The lagrangian density 
ip{X,Y = 0) = X3/2. 



2.3. Extension to non-abelian fields 

The above results can be extended to generalized non-abelian gauge groups. Let us summarize 
the main ideas. The whole analysis will appear in a forthcoming publication. In this case the field 
strength components in the Lie algebra are defined as Fa^u = d^Aau — duAa/j, — gJ2bc CabcAbi_iAcu, 
where 1 < a < {N the dimension of the group) and Cabc are the structure constants. The 
lagrangian density is again defined as ipiy~-A{X,Y) but now the field invariants are defined as 
(we are using the ordinary prescription in the calculation of the traces over the gauge group 
indices) 



^ = "2 J2iFa,.Fn = E(^a - Hi) ■,Y = -- Y.{Fa,.FD = 2J2{E, ■ Ha). (23) 

a a a a 

The associated field equations in presence of gravitation take the form 

E VacM yxFr + V>YFr] = 0, (24) 

c 

where Vac^i = ^ac^ fj. — 9 Z^fec CabcAbfj. is the gauge-covariant derivative. Concerning the 
admissibility conditions, it can be shown that they take in the present case the same form as 
those obtained for NED models. For ESS solutions in this non-abelian case {Ea / 0, Ha = 0, Va), 
the A'^ field equations arising from the v = Q component of the field equations (|24p . in the metric 
([3]), can be combined to give 

r2^/Xv9x|y=o = Q, (25) 

where Q = \jYlia Qa is the "mean-square charge", constructed with the integration constants 
Qa arising from each one of the 1/ = components of the field equations and identified as "source 



color charges" . As a consequence, there is a relation one-to-one between the ESS solutions of the 
generalized non-abelian gauge field theory and those of the NED models with the identification 
ipN-AiX,Y = 0) = (Pned{X,Y = 0) and q ^ Q. The associated energy e{Q) is the same as 
in the abelian case, but now there is a degeneration of energy in the sphere of radius Q. As 
a consequence of these equivalences, the integration of the metric from the Einstein equations, 
and the associated gravitational configurations for generalized non-abelian gauge groups can be 
shown to be the same as for the NED models analyzed above, once the previous identifications of 
lagrangian density functions and charges between abelian and non-abelian cases are performed. 

3. Conclusions and perspectives 

The gravitational configurations associated to general NED models minimally coupled to 
Einstein gravity can be qualitatively characterized without explicitly fixing the form of the 
lagrangian density function. In such a characterization the admissibility conditions play an 
essential role. A natural classification of the ESS solutions into several families is performed 
according to their asymptotic and central- field behaviours, which determine the character (finite 
or divergent) of the flat-space energy. Such NED models have been shown to lead to gravitating 
structures which are asymptotically Schwarzschild-like, or with an anomalous asymptotic non 
Schwarzschild-like behaviour (although still asymptotically flat and well behaved). In both cases 
the structure of horizons was studied, and explicit examples of lagrangian densities belonging 
to the different families of models were given. These procedures are not restricted to the 
asymptotically flat solutions of the Einstein gravity, but can be extended to (Anti-)de Sitter 
spaces and to higher-order curvature gravity theories (e.g. Gauss-Bonnet), a work in progress. 
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